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ABSTRACT 

We use a relativistic ray-tracing code to analyze the X-ray emission from a pressure-supported oscillating 
relativistic torus around a black hole. We show that a strong correlation exists between the intrinsic frequencies 
of the torus normal modes and the extrinsic frequencies seen in the observed light curve power spectrum. This 
correlation demonstrates the feasibility of the oscillating-torus model to explain the multiple peaks seen in black 
hole high-frequency quasi-periodic oscillations. Using an optically thin, monochromatic emission model, we 
also determine how a relativistically broadened emission line and the amplitude of the X-ray modulations are 
dependent on the observer's inclination angle and on the torus oscillation amplitudes. Observations of these 
features can provide important information about the torus as well as the black hole. 

Subject headings: black hole physics - accretion disks - X-rays: binaries 
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consistent relativistic-hydrodynamics simulations. 

This Letter is organized as follows: in Section 2 we sum- 
marize the basic dynamical features of the perturbed relativistic 
tori and describe how we apply the classical radiative trans- 
fer equation to a general relativistic accretion model. In Sec- 
tion 3 we present simulated images of the torus along with in- 
stantaneous line profiles and integrated X-ray light curves for a 
range of inclinations and perturbation amplitudes, representing 
the torus initial conditions. We conclude in Section 4 with a 
discussion of the major results and a look towards future work. 

2. TORUS DYNAMICS AND RADIATIVE TRANSFER 

To briefly summarize the basic properties of the oscillating- 
torus model, we recall that we are considering a non-self- 
gravitating perfect-fluid torus orbiting a Schwarzschild black 
hole (Font & Daigne 2002; Zanotti, Rezzolla, & Font 2003). 
The fluid is assumed to be in circular non-geodesic motion and 
the conditions of hydrostatic equilibrium and of azimuthal sym- 
metry allow the relativistic hydrodynamics equations to be re- 
duced to Bernoulli-type equations. These have particularly sim- 
ple solutions when the fluid is assumed to have a constant spe- 
cific angular momentum and if a polytropic equation of state 
is adopted. In this case, the equations of hydrostatic equilib- 
rium can be integrated analytically to yield the rest-mass den- 
sity distribution inside the torus, with the isobaric surfaces co- 
inciding with the equipotential ones. Hereafter we will consider 
tori with constant specific angular momentum and bear in mind 
that more complex distributions introduce only small quantita- 
tive differences [see Montero et al. (2004); Zanotti et al. (2005) 
for details]. 

Once a stationary equilibrium configuration is constructed, 
it is perturbed with the introduction of a small radial velocity 
expressed in terms of the radial velocity for a relativistic spher- 
ically symmetric accretion flow onto a Schwarzschild black 
hole, i.e. the Michel solution (Michel 1972). More specifically, 
we set the initial radial (covariant) component of the fluid 3- 
velocity as v r — r)(v r ) Michcl , and then use the dimensionless 
coefficient t] to tune the strength of the perturbation, obtaining 
an essentially linear response for 77 < 0.06 (Zanotti, Rezzolla, 
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1. INTRODUCTION 

Recent observations with the Rossi X-ray Timing Explorer 
(RXTE) have revealed the existence of high-frequency quasi- 
periodic oscillations (QPOs) in a number of accreting black 
hole binary systems (Strohmayer 2001; McClintock & Remil- 
lard 2005). In an increasing number of these systems, the 
QPOs appear with integer commensurabilities, generally a 2 : 
3 frequency ratio (Miller et al. 2001; Remillard et al. 2002; 
Homan et al. 2005). Since these modulations are expected to 
originate very close to the black hole, they could be used to 
test gravity in strong-field regimes or extract information on the 
black hole properties. 

Over the years, a large number of theoretical models have 
been developed to explain these observations. Some of the 
more popular models explain the QPOs through magnetic 
flares (Galeev, Rosner, & Vaiana 1979; Poutanen & Fabian 
1999), fluid oscillations in thin disks (Okazaki, Kato, & Fukue 
1987; Nowak et al. 1997), geodesic resonances (Stella & Vi- 
etri 1999; Abramowicz & Kluzniak 2001), and trapped fluid 
oscillations in tori around black holes (Rezzolla et al. 2003; 
Rezzolla, Yoshida, & Zanotti 2003; Montero et al. 2004; Lee, 
Abramowicz, & Kluzniak 2004; Zanotti et al. 2005). 

To evaluate the relative strengths and weaknesses of any of 
these models, it is essential to compare directly the predictions 
of the theoretical models with the observations. To this end, we 
have applied the ray-tracing methods described in Schnittman 
& Bertschinger (2004) to a numerical calculation of the non- 
linear dynamics of a relativistic axisymmetric torus described 
in Zanotti, Rezzolla, & Font (2003). With this combined ap- 
proach we have calculated the profile of a relativistically broad- 
ened emission line, the X-ray light curves, and the power spec- 
tra, all as measured by a distant observer. We have demon- 
strated that a strong correlation exists between the intrinsic 
normal-mode oscillations of a pressure-supported torus and the 
extrinsic observables of the X-ray light curves and power spec- 
tra. It should be noted that this approach is rather different from 
the one recently presented by Bursa et al. (2004), where the 
torus was modeled analytically and thus not the result of self- 
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& Font 2003). As noted in Zanotti et al. (2005), the response 
of the torus is largely independent of the type of perturbation 
and different choices lead to the excitation of the same modes. 

With these initial conditions, the equations of relativistic 
hydrodynamics in a Schwarzschild black hole space-time are 
solved using the axisymmetric, general relativistic code de- 
scribed in Zanotti, Rezzolla, & Font (2003) and Zanotti et 
al. (2005). This makes use of a first-order, flux-conservative 
formulation of the equations, which are solved using a high- 
resolution shock-capturing scheme based on an approximate 
Riemann solver. Second-order accuracy in both space and 
time is achieved by adopting a piecewise-linear cell reconstruc- 
tion procedure and a second-order, conservative Runge-Kutta 
scheme, respectively. As the numerical evolution proceeds, the 
density and pressure in the fluid's local rest-frame, as well as 
the coordinate 4-velocity, are tabulated and stored at each point 
in space-time. 

We recall that the introduction of the perturbations trig- 
gers harmonic oscillations of the torus having centrifugal and 
pressure-gradients as the restoring forces. A careful inves- 
tigation of these oscillations has also revealed that there are 
multiple peaks in the power spectrum with frequencies in a 
sequence of integers: 2:3:4 :... (Zanotti, Rezzolla, & Font 
2003; Zanotti et al. 2005). Subsequent perturbative analyses 
have also shown that these oscillations are indeed p modes, be- 
having as trapped waves within the cavity produced by the torus 
and hence having eigenfrequencies in a sequence of small inte- 
gers (Rezzolla, Yoshida, & Zanotti 2003; Montero et al. 2004). 
The striking analogy between the harmonic relation among the 
p-mode eigenfrequencies and the QPOs observed in black-hole 
systems has then led to the suggestion that QPOs could result 
from basic fluid oscillations of a small accretion torus close to 
the black hole (Rezzolla et al. 2003). The existence of such tori 
appears to be a robust feature of global magnetohydrodynamic 
simulations (DeVilliers, Hawley, & Krolik 2003). 

While attractive for its simplicity and for being based on 
global modes of oscillations that are expected to be present in 
realistic accretion discs, the relativistic torus model has so far 
only suggested a property of an orbiting, pressure-supported 
gas itself. It is not intuitively obvious that "intrinsic" modula- 
tions in the fluid hydrodynamics will produce similar "extrin- 
sic" modulations in the observed light curve. Furthermore, it 
is not clear a priori what might be the relationship between the 
phases and amplitudes of the intrinsic and extrinsic oscillations, 
and how these might compare with the X-ray data. 

To explore this relationship and demonstrate how an oscillat- 
ing accretion model produces a corresponding oscillating signa- 
ture in the observed X-ray light curve, we have calculated the 
trajectories of photons from a distant observer to the emission 
region around a black hole following the methods described 
in Schnittman & Bertschinger (2004). In this approach, the 
photon positions and momenta are tabulated along each ray's 
path in order to recreate a simulated, time-varying image of the 
accreting gas. Given the photon's 4-position and 4-momentum 
along the entire path, the observed spectrum is calculated for 
that ray by integrating the radiative transfer equation with a 
model for the gas emissivity and absorption. 

More precisely, we start from the classical radiative transfer 
equation (Rybicki & Lightman 1979) 

—r- = iu - ol v 1 v , (1) 
as 

where ds is the differential path length and I v , j v , and a v are re- 
spectively the spectral intensity, emissivity, and absorption co- 



efficient of the fluid at a frequency v. These variables are typi- 
cally defined in the rest-frame of the gas as functions of its local 
temperature and density. We have explored a number of differ- 
ent emission models, which will be presented in greater detail 
in a companion paper. For simplicity, the results discussed here 
are primarily based on an optically thin gas emitting isotrop- 
ically and monochromatically at frequency z/ em with a v = 
and j v oc p5(v - z/ em ). 

We then incorporate relativistic effects by defining a local 
orthonormal tetrad at each point along the integration path (for 
the Schwarzschild geometry, this is simply the coordinate sta- 
tionary observer frame). Once transformed to this locally-flat 
tetrad, only special relativistic effects must be included. Fol- 
lowing Rybicki & Lightman (1979), equation (1) becomes 

where primed and unprimed variables are measured in the rest- 
frame of the gas and in the stationary tetrad, respectively. The 
special relativistic Doppler shift between the photon path and 
the fluid can be written as 

^= 7 (l-/?cos0), (3) 

where (3 = v/c, 7 = 1/ \/l — (3 2 , and 9 is the angle between 
the photon direction and the gas, all measured in the stationary 
tetrad. In addition, as the photon bundle propagates through the 
global curvature around the black hole, the spectral intensity 
at a given frequency evolves as the photons are gravitationally 
red-shifted, maintaining the Lorentz invariance of I v jv 3 . 

3. LIGHT CURVES AND POWER SPECTRA 
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FIG. 1. — Normalized X-ray light curves from oscillating tori with incli- 
nations i = 10° (solid) and 70° (dashed), for a perturbation amplitude of 
r\ = 0.02. The local minima in the light curves correspond to minima in the 
torus size, when it is closer to the black hole. The inset shows the same light 
curves in greater detail. 

Figure 1 shows the light curves for a torus orbiting around a 
Schwarzschild black hole of mass M = 10 M Q at inclination 
angles of 10° and 70°. The amplitude of the initial perturbation 
was rj — 0.02, but qualitatively similar behaviors are seen for 



QPOs in the Light Curves from Oscillating Tori 



3 




E obs /E em E o6 /E 1In E „b,/ E ,m 

FIG. 2. — Upper: Ray-traced images of the oscillating torus, shown at various phases of the fundamental p mode, with a logarithmic color scale of the X-ray 
intensity. Lower: Broadened emission line spectrum for each frame. The triple-peak spectrum is caused by relativistic beaming towards and away from the observer, 
along with a central peak due to gravitational lensing of the far side of the accretion torus. The inclination is i = 70° and the perturbation amplitude is r] = 0.06. 



perturbations in a linear regime, i.e. 81 /I oc 77 for 77 < 0.06, 
where I = j Q I v dv. Although it is clear that the light curve 
has a quasi-periodic behavior and we find that this is strictly 
related to the oscillating behavior of the torus at the same fre- 
quencies, the modulation of the intensity is the combined re- 
sult of several different relativistic effects. In particular, for 
the optically thin emission model, the minimum of the inten- 
sity is reached in the "compression" phase of the oscillation, 
when the torus size is smaller and thus is closer to the black 
hole. In this case, the gravitational red-shift reduces the ob- 
served photon energies (as well as photon number through the 
invariance of I„/v 3 ) as the light has to escape from a deeper 
potential. At the same time, the intensity is also varied by 
the special relativistic beaming of photons emitted towards and 
away from the observer. Finally, smaller contributions to the 
intensity modulation also come from the transverse ("second- 
order") Doppler shift, and from the gravitational lensing of the 
far side of the torus, magnifying a small region of emission that 
is moving transverse to the observer (Beckwith & Done 2004; 
Schnittman, Homan, & Miller 2005). 

Because of the conservation of angular momentum during 
the oscillations, the fluid velocities increase as the torus ap- 
proaches the black hole, thus enhancing the relativistic beaming 
of photons toward the observer. For large inclination angles (i.e. 
when the torus is almost "edge-on"), this beaming is particu- 
larly intense and serves to compensate for the intensity decrease 
due to the gravitational red-shift. As a result, the two major rel- 
ativistic effects counter each other and the intensity modulation 
is smaller. For small inclination angles (i.e. when the torus is 
almost "face-on"), the beaming and gravitational lensing can 
largely be ignored and the intensity modulation is dominated 
by the gravitational red-shift and is thus relatively larger (com- 
pare dashed and solid lines in Fig. 1). For these smaller incli- 
nation angles, simple estimates can be made through the modu- 
lated red-shift (5^ b s = t'obs — (^obs) an d from the invariance of 
Iv/i/ 3 ,to obtain 5I/I w 3(<5^ b s /t'obs)- The red-shift variation 
^obs/^obs is in turn linearly proportional to 77. 

It is important to underline that the dependence of the mod- 
ulation of the light curve on the inclination angle is the qualita- 
tive opposite of the hot spot model described in Schnittman & 



Bertschinger (2004) and could serve to distinguish between the 
two models as more observations become available. However, 
because this dependence is a function of the emission model 
used, further studies are necessary. 

In the upper row of Figure 2 we show three snapshots of an 
oscillating torus at different phases of a single period P, with 
t/P = 0.0, 0.25, 0.5, for a torus with inclination angle i = 70° 
and a perturbation amplitude 77 = 0.06. In the lower row, we 
show a broadened emission line spectrum for each frame with 
the typical "multi-horned" features, where the two main peaks 
are caused by photons emitted towards and away from the ob- 
server and thus being blue- and red-shifted respectively. The 
smaller intermediate peak, on the other hand, is due to the grav- 
itational lensing of the far side of the torus (Beckwith & Done 
2004; Schnittman, Homan, & Miller 2005). 

When the torus is more compact, closer to the black hole, 
and the fluid velocities are comparatively larger, the observed 
emission line is widest as the red-shifted wing moves to slightly 
lower energies and the blue-shifted wing to slightly higher ones 
(see panel at t/P = 0.0 in Fig. 2). Furthermore, because of 
the relativistic invariance of I„/^ 3 , there are more blue-shifted 
photons, which increases the total observed flux, compensating 
in part for the gravitational red-shift of the smaller torus. On 
the other hand, when the torus is at its maximum size and far- 
ther away from the black hole, the blue- and red-shifts in the 
line spectrum are smaller, but the height of the peaks is com- 
paratively larger (see panel at t/P = 0.5 in Fig. 2). 

To quantify more precisely the "quasi-periodicity" in the 
light curve expected from an oscillating torus, we show in Fig- 
ure 3 the power spectra from the light curves reported in Figure 
1, in units of [(rms/mean) 2 Hz -1 ]. The power is clearly dom- 
inated by a peak at ~ 58 Hz, which coincides with the lowest 
p mode of the torus (indicated as the fundamental /). Further- 
more, the power spectra also show smaller peaks at the over- 
tones o\ ~ 3/2/, 02 ~ 2/ and 03 ~ 5/2/, thus demonstrating 
that the light curve possesses the same harmonic behavior of 
the underlining hydrodynamics. 
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FIG. 3. — Power spectra of the light curves plotted in Figure 1. Note the 
harmonic ratios among peaks, with o\ ~ 3/2 /, 02 ~ 2 /, and 03 ~ 5/2 /. 

To distinguish between transient modes and the more per- 
sistent fundamental oscillations, we have also calculated power 
spectra for longer time series, corresponding to more than 30 
periods of oscillations. In this case we find that the power in 
the higher modes is reduced and may even disappear if the time 
series is restricted to the later stages of the oscillations. This 
result, which is a numerical artifact for low-amplitude pertur- 
bations, could however serve as a guide in interpreting the ob- 
servations. In a real accreting system, in fact, it is likely that 
any perturbation lasts only a short time, as the turbulent condi- 
tions of the accretion disk would continually create and destroy 
coherent oscillating tori. As a result, these short-lived pertur- 
bations could maintain more significant power in the higher 
harmonics than what was produced with these simulations. In- 
deed, a characteristic torus lifetime of only 3 — 5 periods is 
sufficient to explain the observed oscillator quality factors of 
Q w 10 — 15 (Schnittman 2005). Furthermore, the discontinu- 
ous phase shift between subsequent tori provides then a natural 
explanation for the broadening of the QPO peaks and the small 
changes in frequencies observed as a "jitter" in the peaks. 

4. DISCUSSION AND CONCLUSIONS 

We have demonstrated the positive correlation between the 
intrinsic normal mode oscillations of a pressure-supported torus 
and the extrinsic X-ray light curves and power spectra as seen 



by a distant observer. This confirms the feasibility of the 
oscillating-torus model as an explanation for the integer ratios 
seen in high-frequency QPO peaks. The specific parameters of 
the torus model still require further investigation in order to best 
fit the QPO data. 

For the simple emission model considered here, the varia- 
tion in the light curve is primarily caused by the gravitational 
red-shift of photons coming from different radii as the torus 
moves in and out of the black hole's potential well. Unlike 
the relativistic hot spot model, the oscillating torus model pre- 
dicts higher amplitude variations in the light curve for smaller 
inclination angles, while at higher angles the special relativis- 
tic beaming and gravitational lensing counter the gravitational 
red-shift, reducing the variations in flux. This difference could 
be crucial to distinguish the two models as more observations 
become available. 

The results reported in this paper serve essentially to demon- 
strate an issue of principle and have therefore been restricted to 
a simplified scenario. Even with this advanced analysis tool, 
we are still limited by the inherent uncertainty in the emis- 
sion mechanisms and the geometry of the surrounding accretion 
disk and corona. The fact that the QPOs tend to appear exclu- 
sively in the Steep Power Law ("Very High") state of the black 
hole and are most significant in the higher energy channels on 
RXTE (McClintock & Remillard 2005) suggests that the spe- 
cific emission model may be constrained by the features of the 
X-ray energy spectrum. These spectral features are also impor- 
tant in characterizing the properties of the surrounding thermal 
disk and hot corona. 

Work is in progress to make the results presented here more 
realistic by analyzing the dynamics of the torus in a Kerr space- 
time, by considering electron scattering through a corona, and 
by including optically thick line emission and thermal free-free 
emission with a Kramer's opacity law. Initial calculations show 
qualitatively similar results for the line emission model, but the 
thermal emission model predicts a "phase inversion" with re- 
spect to the light curves shown in this paper. This is due to the 
fact that for a barotropic equation of state the temperature in- 
creases at smaller radii when the density and pressure increase. 
This higher temperature produces greater X-ray emission, out- 
weighing the gravitational red-shift of the smaller radius. These 
results will be presented in greater detail in a forthcoming pa- 
per. 

JDS is grateful to Edmund Bertschinger for his insights and en- 
couragement. Helpful discussions with Olindo Zanotti are also 
acknowledged. Support comes by NASA grant NAG5- 13306. 
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